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Abstract. We constructed diagnostic diagrams using 
emission line ratios and equivalent widths observed in 
several independent samples of H ii galaxies. Significant 
trends are seen, both in the line ratio diagrams, and in di- 
agrams relating line ratios to the equivalent width of H/3. 
The diagrams are compared to predictions from photoion- 
ization models for evolving starbursts. This study extends 
the work of Stasinska & Leitherer (1996) by including ob- 
jects with no direct determination of the metallicities, and 
by using updated synthesis models with more recent stel- 
lar tracks and atmospheres. 

We find that H ii galaxies from objective-prism sur- 
veys are not satisfactorily reproduced by simple models 
of instantaneous starbursts surrounded by constant den- 
sity, ionization bounded H ii regions. The observed re- 
lations between emission line ratios and H/3 equivalent 
width can be understood if older stellar populations gen- 
erally contribute to the observed optical continuum in H ii 
galaxies. In addition, different dust obscuration for stars 
and gas and leakage of Lyman continuum photons from 
the observed H ii regions can be important. As a result, 
H II galaxies selected from objective-prism surveys are not 
likely to contain significant numbers of objects in which 
the most recent starburst is older than about 5 Myr. This 
explains the success of the strong line method to derive 
oxygen abundances, at least in metal poor H ii galaxies. 

The observed increase of [O i] /H/3 with decreasing H/3 
equivalent width can result from the dynamical effects of 
winds and supernovae. This interpretation provides at the 
same time a natural explanation of the small range of 
ionization parameters in giant H ii regions. The classi- 
cal diagnostic diagram [O ill] /H/3 vs [O ii]/H/3 cannot be 
fully understood in terms of pure photoionization models. 
The largest observed [O ii]/H/3 ratios require additional 
heating. 

The [N ii]/[0 ii] ratio is shown to increase as the H/3 
equivalent width decreases. A possible explanation is an 
N/0 increase due to gradual enrichment by winds from 
Wolf-Rayet stars on a time scale of ~ 5 Myr. Alternatively, 
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the relation between N/0 and 0/H could be steeper than 
N/0 oc 0/H°-^, with a previous stellar generation more 
important at higher metallicities. 

Key words: Galaxies: abundances - Galaxies: evolution 
- Galaxies: ISM - Galaxies: starburst - Galaxies: stellar 
content - ISM: H ii regions - Stars: Wolf-Rayet 



1. Introduction 

Isolated extragalactic H ii regions, which we will simply 
refer to as H ii galaxies, are powered by clusters of hot, 
massive stars that ionize their environment. Recent analy- 
ses of such objects, comparing the observed nebular emis- 
sion lines, colors, and stellar features with evolutionary 
models (Mas-Hesse & Kunth 1991, 1999, Schaerer et al. 
1999, Stasinska & Schaerer 1999) found that most have 
experienced a recent, quasi instantaneous burst of star for- 
mation. In addition, there is growing evidence for an older 
stellar population (Telles & Terlevich 1997, Papaderos et 
al. 1998, Raimann et al. 2000a). Some clues regarding the 
evolution of giant H ii regions, their properties and stel- 
lar content may be obtained by considering H ii galaxies 
in different evolutionary stages with respect to the most 
recent star formation event. 

The present paper follows this approach: we compare 
the emission line properties of a large sample of H II 
galaxies with those of photoionization models in which 
the ionizing continuum is provided by synthetic models of 
evolving starbursts. In a previous paper (Stasinska & Lei- 
therer 1996, hereinafter SL96), we considered a sample of 
metal poor galaxies in which the oxygen abundance could 
be derived directly from observations using the electron 
temperature sensitive [O ill] A4363 line. Thus the derived 
abundances were unambiguous and model independent. 
Unfortunately, we were automatically restricted to young 
starbursts by this requirement. As the most massive stars 
disappear, the radiation field gradually softens inducing a 
decrease of the [O III] emission, and the weak [O III] A4363 
becomes undetectable after about 5 Myr. Alternatively, 
the oxygen abundances can be estimated via strong line 
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methods (Pagcl ct al. 1979, McGaugh 1991, 1994). The 
relevance and accuracy of these methods is however still 
under debate. We address this question in Sect. 4.2. 

Lifting the restriction to the youngest starbursts im- 
mediately leads to inclusion of objects with a priori un- 
known metallicities. However, the distribution in metallic- 
itics should be the same for old and young objects, per- 
mitting an interpretation of diagnostic diagrams based on 
emission lines. This advantage of a large sample of bona 
fide young and old H ii galaxies is explored in the present 
work. With respect to SL96, we also use updated evolu- 
tionary synthesis models with more recent stellar tracks 
and stellar atmospheres as input for our photoionization 
models. 

The study in the present paper is complementary to 

work on giant H ii regions in spiral galaxies, such as pub- 
lished by Garcia- Vargas & Diaz (1994), Garcia- Vargas et 
al. (1995), and more recently Bresolin et al. (1999) and 
Dopita ct al. (2000). First, our observational samples arc 
based on extragalactic H ii regions, for which the equiva- 
lent width of H/3 can be used as a direct first-order indica- 
tor of the age of the latest starburst (Dottori 1981). Sec- 
ond, our samples are biased towards dwarf galaxies with 
correspondingly low average metallicities. 

In Section 2 we present the observational samples, dis- 
cuss possible selection effects, and display a series of ob- 
servational diagrams. In particular, we show that emission 
line trends arc seen not only in classical emission line ratio 
diagrams, but also in diagrams relating emission line ra- 
tios with the H/3 equivalent width. In Section 3 we describe 
the models, present our model grid, and show theoretical 
diagrams for evolving starbursts. In Section 4 observations 
and models are compared. The main conclusions are sum- 
marized in Section 5. 

2. Sample selection 

2.1. The database 

The prerequisites for our observational data base are the 
following. We require a large sample of objects with uni- 
form observation parameters. The sample must be large 
since there arc several factors determining the line spec- 
tra, the most important being metallicity and starburst 
age. We wish to use information from the He i A5876 line, 
which is directly related to the mean effective tempera- 
ture of the ionizing stars (cf. Stasinska 1996). Deep spec- 
troscopy is required for this generally weak line. In ad- 
dition, we want to make use of a direct estimator of the 
starburst age, such as provided by the equivalent width of 
H/3 (hereafter EW{ILP)). While this estimator has its well- 
known drawbacks (effects of dust, incomplete absorption 
of the stellar ionizing photons by the observed nebula), it 
is relatively independent of the assumed nebular proper- 
ties, as compared with emission line ratios such as [S ii] 
AA6717-h30/H/3or [S ii] AA6717-F30/[S iii] A9069 proposed 



by Garcia- Vargas et al. (1995). However, EW(HP) is diffi- 
cult to use for giant H ii regions with a strong underlying 
old stellar component, such as in nuclear starbursts or gi- 
ant H II regions in evolved galaxies. Our sample is thus 
restricted to giant H ii regions in irregular or blue com- 
pact galaxies, where the old stellar component is expected 
to be weaker. 

The largest, homogeneous sample meeting these re- 
quirements is the spectrophotometric catalogue of H ii 
galaxies by Terlevich et al. (1991). It contains spectra of 
425 emission line galaxies, among which about 80% are 
classified as genuine H ii galaxies by the authors (as op- 
posed to active galactic nuclei and nuclear starbursts). 
Further imaging of these objects identified some of these 
H II galaxies as giant H ii regions in spiral galaxies (Telles, 
private communication). Those objects were removed to 
produce the final sample, hereafter referred to as the Ter- 
levich sample, which consists of 305 objects. Most of the 
observations for the Terlevich catalogue were made with 
2-m class telescopes, and the typical signal-to-noise in the 
continuum is 5. Note, however, that second-order contam- 
ination could affect some of the line ratios for A > 6000 A 
(Terlevich et al. 1991). 

As a control sample, we consider a much more re- 
stricted collection of isolated extragalactic H ii regions 
with very high quality observational data. These are the 
blue compact galaxies observed by Izotov and coworkers 
(Izotov et al. 1994, 1997, Thuan et al. 1995, Izotov & 
Thuan 1998, Guseva et al. 2000) on 2-m and 4-m class 
telescopes. Their signal-to-noise in the continuum is typi- 
cally 20 - 40. The latter sample, referred to as the Izotov 
sample, comprises 69 objects. 

Recently, another survey of emission line galaxies 
meeting our requirements has been published. This is the 
catalogue of Popescu & Hopp (2000) which contains 90 ob- 
jects, out of which 70 are classified as H ii galaxies. The 
quality of the data is not as good as for the Izotov sample, 
but the sample is more homogeneous, as explained below. 

2.2. Selection effects 

Objective-prism surveys of emission line objects tend to 

underestimate the true proportion of objects with weak 
emission line equivalent widths. On the other hand, gen- 
uine, isolated giant extragalactic H ii regions will always 
be discovered in such surveys if they are bright enough. 
The Terlevich catalogue comprises 75% (25%) of the non- 
quasar emission line objects identified in the Michigan 
(Tololo) objective-prism survey, and additional biases are 
not expected to be important. Therefore, one expects a 
roughly uniform distribution of ages of the most recent 
starbursts, and the same distribution of metallicities at 
each age bin above a certain threshold in H/3 equivalent 
width in the Terlevich sample. 

In the case of the Izotov sample, the objects are mainly 
blue compact galaxies from the First and Second Byu- 
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rakan objective-prism surveys. The sample is obtained 
by merging two sub-samples. One is composed of ob- 
jects which had been selected for their low metallicities 
{Z < Z0/IO) as estimated from low signal-to-noise 6-m 
telescope spectra. The other subsample (the 39 objects 
studied by Guseva et al. 2000) is composed of objects se- 
lected on the basis of broad Wolf-Rayet features seen in 
the 6-m telescope spectra. The metallicities in the latter 
subsample range between Zq/AO and ZqQ The H/3 equiva- 
lent widths in the Guseva subsample are, on average, lower 
than those in the other Izotov et al. subsample, and con- 
tamination from old stellar populations is likely present at 
least for some objects known as nuclear starbursts. Com- 
pared to the Terlevich sample, the Izotov sample is af- 
fected by strong biases which must be kept in mind when 
interpreting the diagnostic diagrams shown below. 

The Popescu & Hopp sample is a complete subsam- 
ple of emission line galaxies found on the objective-prism 
plates for the Hamburg Quasar Survey. The single selec- 
tion criterion was the location of the galaxies with respect 
to voids. Therefore, as for the Terlevich sample, we ex- 
pect a roughly uniform distribution of ages, and the same 
metallicity distribution at each age bin. The sample is 
however much smaller than the Terlevich sample; there- 
fore statistical fluctuations may not be negligible. 

2.3. Reddening corrections 

Izotov and coworkers used an iterative procedure to dered- 
den their sample where the reddening constant and stellar 
absorption at H/3 are determined simultaneously. Such a 
procedure was feasible due to the high signal-to-noise of 
the data. They assumed the Whitford (1958) extinction 
law and took for the intrinsic IIa/H/3 ratio the recom- 
bination value corresponding to the electron temperature 
derived from [O ill] 4363/5007, when available. In most 
cases, the contribution of the stellar absorption in H/3 
turns out to be only a few percent, except in some ob- 
jects from the Guseva sample which includes a few nuclear 
starbursts (see Guseva et al. 2000, Schaerer et al. 2000). 

The line intensities in the Popescu & Hopp (2000) sam- 
ple are published corrected for reddening. For the galax- 
ies with a strong underlying continuum relative to the H/3 
emission {EW(liP) < 20 A), they assumed a constant 
underlying absorption at H/3 of 2 A. The reddening was 
computed from Ha/H/3 assuming an intrinsic ratio of 2.87 
and using the Howarth (1983) extinction law. 

In the case of the Terlevich sample, the data are pub- 
lished without reddening correction. We have dereddened 
them adopting the Seaton (1979) reddening law and an 
intrinsic Hq;/H/3 ratio of 2.87. We assumed a constant un- 
derlying stellar absorption of 3 A for all the spectra in the 

^ For 9 objects the electron temperature was too low for tem- 
perature sensitive ratios to be measured, and the metallicity 
was estimated with the strong line method (Pagel et al. 1979) . 



Ha and H/3 lines, which is predicted by models for syn- 
thetic absorption line spectra in starburst galaxies (Olofs- 
son 1995a, Gonzalez Delgado et al. 2000). Unfortunately, 
spectra in the Ha region have second order contamination 
in many objects (cf. Terlevich et al. 1991), affecting the 
reddening correction. Alternatively, one could use H7/H/3 
instead of Ha /H/3 to determine the reddening correction, 
but this procedure is plagued by the fact that H7 is a much 
weaker line than Ha, introducing substantial errors due to 
the reduced signal-to-noise, and by the lack of published 
H7 equivalent widths. In any event, for objects without 
published Ha data, we used H7/H/3 to derive the redden- 
ing, assuming that the stellar (negative) contribution to 
the observed H7 line is three times as large as for the H/3 
line. 

The reddening correction introduces some uncertain- 
ties into each sample. Even in the case of the Izotov sam- 
ple, where dereddening was performed in a completely 
self-consistent way, the reddening correction may not be 
perfect. First, departures from the standard reddening law 
exist (e.g. Mathis & CardeUi 1988 or Stasinska et al. 1992). 
Second, the intrinsic Ha/H/3 ratio in a nebula is not neces- 
sarily the recombination value. In objects with high elec- 
tron temperatures, collisional excitation increases Ha/H/3 
with respect to the recombination value, in a proportion 
that depends on the quantity of neutral hydrogen in the 
emitting regions (Davidson & Kinman 1985, Stasinska & 
Schaerer 1999). Adopting the same underlying absorption 
of 3 or 2 A at H/3 for all the objects of the Terlevich and 
Popescu & Hopp samples may not be justified if stellar 
populations from previous star-forming events contribute 
significantly to the continuum (see below). Among the line 
ratios discussed in this paper, the one most affected by 
reddening is [O 11] /H/3. It is prudent to assume that the 
[O 11] /H/3 ratios may well be uncertain by up to about 
30%, and even more for some objects in the Terlevich sam- 
ple. 

2.4. Observed emission line trends 

The nine panels of Figs. 1, 2, and 3 are various obser- 
vational diagrams for the Terlevich, Izotov, and Popescu 
& Hopp samples, respectively. Panels a through f in these 
figures show the behavior of various emission line ratios as 
a function of i?VF(H/3). In spite of some (expected) disper- 
sion, striking line trends are seen. These trends are very 
similar in the three samples, and are shown for the first 
time. 

Some differences among the samples are likely due to 
the quality of the data. The smaller dispersion in He i 
A5876/H/3 in the Izotov sample, for example, is probably 
due to the better signal-to-noise. We also note that the 
Terlevich and Popescu & Hopp samples contain some ob- 
jects with rather large [O li]/H/3: 20 % of the Popescu 
& Hopp sample and 3% of the Terlevich sample have 
[O ii]/H/3 >5, while no object of the Izotov sample has 
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such large values for this ratio. However, if we retain only 
the 50 objects from the Terlevich sample with measure- 
ments in the four lines [O ii] A3727, [O III] A5007, He i 
A5876 and [N ii] A6584, the diagrams become very similar 
to those corresponding to the Izotov sample. Therefore, we 
will ignore the objects with [O ii]/H/3 > 5 since such high 
ratios may be attributed to poor signal-to- noise and/or to 
inadequate reddening corrections. 

We generally note a gradual increase in the disper- 
sion and a decrease of the average value of [O iii]/H/3 as 
EW{B.P) decreases, while the He i A5876/H/3 ratio re- 
mains remarkably constant. The [O l]/H/3 ratio, on the 
other hand, steadily increases with decreasing EW{HP), 
extending the trend already noted by SL96 to lower 
EWinp). The [O ii]/H/3 increases and tends to level off 
at EWinp) about 30 A. The ([O n]+ [O iii])/H/3 ratio 
is rather constant at the largest EW(H(3), and becomes 
more dispersed with a tendency to decrease as EW(Rf3) 
decreases. The [N ii]/[0 ii] ratio shows a clear tendency 
to increase as EW(lif3) decreases, although the dispersion 
is larger than in the remaining diagrams. Overall, all the 
trends of emission line ratios with H/3 equivalent width 
are significant and require an explanation. 

Finally, panels g, h, and i in Figs. 1-3 show the 
standard emission line ratio diagrams ([O iii]/H/3 vs 
[O ii]/H/3, [O iii]/H/? vs [S ii]/H/3 and [N ii]/[0 ii] vs 
([O ll]-|-[0 iii])/H/3), similar to those that have been 
widely used in the literature to distinguish H ii regions 
from active galaxies (Baldwin et al. 1981, Veilleux & Os- 
terbrock 1987, van Zee et al. 1998, Martin & Friedh 1999) 
and for studies of abundance ratios in H II regions (Mc- 
Gaugh 1994, Ryder 1995, Kennicutt & Garnett 1996, van 
Zee et al. 1998, Bresolin et al. 1999). For comparison, we 
have plotted in Fig. 4 the same three emission line ra- 
tio diagrams for the sample of giant H ii regions in spiral 
galaxies from McCall et al. (1985). 

The H II galaxies of our samples clearly appear like an 
extension of the giant H ii region sequence of McCall et 
al. toward the high [O iii]/H/3 end and towards the low 
[N ii]/[0 ii] end. Some of the H ii galaxies (mainly in the 
Terlevich sample) reach values of [N ii]/[0 ii] which are as 
high as those observed in the giant H ii regions of the cen- 
tral parts of spiral galaxies. The true scatter in [O ii]/H/3 
is probably smaller than it appears in our figures: as dis- 
cussed in Section 2.3, [O ii]/H/3 is strongly affected by 
uncertainties in the dereddening procedure]^ Note the 
pronounced bending of the [O iii]/H/3 vs [O ii]/H/3 and 
[O iii]/H/3 vs [S ii]/H/3 sequences toward the left at high 
[O iii]/H/3 in the three H ii galaxies samples. The relation 
between [N il]/[0 ii] and ([O u] + [0 iii])/H/3 is extremely 
tight in all the samples we consider. This was also seen 
in diagrams originally studied by the authors mentioned 



^ Plotting [O iii]/H/3 vs [O ll]/H/3 from the Terlevich sample 
using line ratios not corrected for reddening results in a smaller 
scatter. 



above. Considering that line ratios in nebulae are deter- 
mined by several factors (abundance ratios, intensity and 
spectral distribution of the ionizing radiation field, density 
distribution of the nebular gas), the question arises as to 
why the observed sequences are so narrow. 

We will try to understand the observational data in 
the light of the photoionization models described in the 
next section. 



3. The models 

3.1. Technique and assumptions 

The evolutionary synthesis models of Schaerer & Vacca 
(1998) are used to predict the theoretical spectral energy 
distributions^. They are based on the non-rotating Geneva 
stellar evolution models, with the high mass-loss tracks of 
Meynet et al. (1994). The spectral energy distributions for 
massive main-sequence stars are those given by the CoStar 
models (Schaerer & de Koter 1997) taking into account 
the effects of stellar winds, non-LTE, and line blanketing. 
The pure He models of Schmutz et al. (1992) are used for 
Wolf-Rayet stars. The spectral energy distributions from 
the plane-parallel LTE models of Kurucz (1991) arc used 
for the remaining stars which contribute to the continuum. 

To illustrate the differences between the ionizing spec- 
tra adopted here with the recent Starhurst99 models (Lei- 
therer et al. 1999) and the earlier Leitherer & Heckman 
(1995) models adopted in SL96, we plot the ratio of the 
He*^ and Hc+ ionizing photons (Qhc, Qhc+ respectively) 
to Qh, the number of H ionizing photons issued from var- 
ious models in Figure || Differences with the Leitherer 
& Heckman (1995) models (left panel) are mostly due 
to the updated stellar tracks, which in particular, lead 
to a shorter Wolf-Rayet (WR) phase in the solar metal- 
licity model shown, and thus to a more rapid decline of 
the hardness of the ionizing spectrum. From this figure 
it is also apparent that at solar (and higher) metallici- 
ties the changes brought about by recent stellar tracks 
have a rather important impact on the ionizing fluxes, 
which affects photoionization models based on the older 
Leitherer & Heckman (1995) spectra (e.g. SL96, Bresolin 
et al. 1999). As both the Schaerer & Vacca (1998) and 
Starburst99 models use the same stellar tracks, the differ- 
ences (right panel) in the predicted spectra are only due 
to different treatments of stellar atmospheres. Outside the 
phases with WR stars (i.e. at ages ^ 3 Myr or ^ 6 Myr), 
the Schaerer & Vacca spectrum in the He" ionizing con- 
tinuum is somewhat harder due to the use of improved O 
star model atmospheres (cf. also Schaerer & Vacca 1998, 
Fig. 5). Differences in the He^" ionizing flux (non- negligible 
during part of the WR phase; cf. Schaerer & Vacca) are 
due to slightly different assumptions on the link between 



The predic tions from these synthesis models are available 
on the Web at http://webast.ast.obs-mip.fr/people/schaerer/. 
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the stellar evolution and atmosphere models and other 
numerical differences. 

As also apparent from Fig. ^, WR stars are predicted 
to provide a non- negligible fraction of the ionizing flux, es- 
pecially at high energies during the WR-rich phase, whose 
duration increases with metallicity (cf. Schaerer & Vacca 

1998) . The predicted ionizing fluxes are likely more uncer- 
tain during these phases, due to 1) the neglect of metals 
in the WR atmosphere models of Schmutz et al. (1992), 
and 2) uncertainties in relating the stellar interior and at- 
mosphere models (cf. review by Schaerer 2000). However, 
while the predicted spectra are possibly too hard at high 
metallicities (see also Crowther et al. 1999, Bresolin et al. 

1999) , there are good indications that at low metallicities 
- of main concern here - the non blanketed WR atmo- 
spheres combined with the latest stellar tracks provide 
a reasonable description of the ionizing spectra (e.g. de 
Mello et al. 1998, Luridiana et al. 1999, Gonzalez Delgado 
& Perez 2000, Schaerer 2000 ). 

We used the code PHOTO to build the photoionization 
models, with exactly the same atomic data as in SL96. 

3.2. The model grid 

We built a wide model grid in order to cover the full pa- 
rameter space found in our observational samples^. 

Most models are computed for instantaneous bursts 
with a Salpeter initial mass function and an upper stel- 
lar limit Mup = 120 Mq. The evolution of the burst is 
followed in time steps of 10^ yr, starting from an age of 
10'* yr. Such a time step may be too long for a detailed 
description of the Wolf-Rayet stage in the evolution of 
the burst at metallicities lower than Zq/5 (see Schaerer & 
Vacca 1998), but is a good compromise to cover the evolu- 
tion until 10 Myr. A second set of models was computed 
for an upper stellar limit of 30 M0 instead of 120 M0. 
Another set of models was built with Mup = 120 M0 and 
a Salpeter IMF but for a star formation extended over a 
period of 10 Myr. 

We considered five metallicities: 2, 1, 0.4, 0.2 and 0.02 
times the solar metallicity. The metallicities of the stars 
and the nebular gas are assumed to be the same. As in 
SL96, the metallicity is defined as the oxygen abundance 
(oxygen is the main gas coolant). 0/H is taken equal to 
8.51 X 10^'' for the solar abundance set. The abundances 
of the other elements relative to oxygen are given by the 
prescription of McGaugh (1991). In particular, Hc/H ~ 
0.0772 + 15 (0/H) and log N/0 = 0.5 log (0/H) + 0.4. 
The ionized gas is assumed to be dust free. 

The third main parameter that determines the emis- 
sion line ratios and equivalent widths of a model is the 
gas distribution. The photoionization models in our grid 



Tables giving infrared, optical and ultraviolet line intensi- 
ties for all t he models discussed in this paper are available on 



are spherically symmetric with constant density, and are 
ionization bounded. Most models presented in this paper 
are spheres uniformly filled with gas at a density of n = 
10 cm~'^. We also constructed some models represented 
by hollow spheres with / = 3 and n = 200 cm~^, where 
/ is the inner radius in units of Stromgren radius for a 
full sphere (see Stasihska & Schaerer 1997). These two ex- 
treme cases are intended to highlight the sensitivity of the 
results on the adopted density distribution. It can easily 
be shown that such hollow sphere models have the same 
ionization parameter (averaged over the nebular volume) 
as the filled sphere models with n = 10 cm^"^ and a total 
ionizing luminosity 1000 times lower. 

Four values of the ionization parameter are consider- 
ered: those that would be produced by a star cluster of 
M*= 1, 10'^, 10^, and 10^ Mq^, providing a range in ion- 
ization parameters of a factor 1000 at a given epoch. One 
should realize that our choice of density distribution is 
quite arbitrary. So far, little is known theoretically or ob- 
servationally about the systematics of the evolution of the 
nebular geometry with time in giant H 11 regions. Our 
model sequences are not necessarily meant to describe the 
true evolution of a giant H 11 region with time. For ex- 
ample, it may well be that, due to the combined effects 
of expansion or stellar winds, the ionization parameter of 
giant H 11 regions decreases with time more than by the 
sole effect of the decrease of the number of ionizing pho- 
tons. This was suggested by SL96 from the comparison of 
a sample of young H 11 galaxies with their grid of models. 
Furthermore, the time variation of the ionization param- 
eter depends on the adopted evolution of the geometry. 
For example, the ionization parameter U in our models 

1 /3 

varies as (Qh being the number of stellar photons 

able to ionize hydrogen per unit time) while in the models 
of Garcia- Vargas & Diaz (1994) and Garcia- Vargas et al. 
(1995), who constructed shells with a fixed outer nebular 
radius, U varies as Qh- 

Figures 6,7 and 8 correspond to our full sphere model 
sequences for instantaneous bursts with Afup — 120 Mq 
and = 10^, 10^, and 10^ Mq, respectively. Figure 9 
corresponds to a Mi, — 10^ Mq burst uniformly extended 
over a period of 10 Myr. (In this case the time step is 2 
Myr instead of 1 Myr.) Figure 10 is for a M^, = 10^ burst 
with Mup = 30 Mq. Figure 11 corresponds to our hollow 
sphere model sequences with Mup = 120 Mq and M^, = 
10^ Mq. Each model is represented by a symbol: circle for 
Z/Zq = 2, square for Z/Zq = 1, triangle for Z/Zq = 0.4, 
diamond for Z/Zq = 0.2 and plus sign for Z/Zq = 0.05, 
and the successive models in a sequence are linked by a 
thin line. 

It is important to realize that, while the model predic- 
tions for optical forbidden lines like [O ill] A5007, [O 11] 
A3727, [N 11] A6584 are relatively robust and accurate for 
metallicities below solar, the situation worsens as one goes 
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to higher nictaUicitics. The reason is that, as metaUicity 
increases, cooling by heavy elements becomes more and 
more important, and is gradually shifted from the optical 
forbidden lines to the fine structure, infrared lines, which 
have lower excitation potential. The emissivities of the in- 
frared lines are little dependent on the electron tempera- 
ture. As a consequence, the computed value of the electron 
temperature, obtained by equating the energy gains and 
losses of the electrons, strongly depends on the computed 
energy loss. Therefore, any small variation in the energy 
loss will induce an appreciable change in the emission of 
the optical forbidden lines. For example, since the infrared 
fine structure lines are easily collisionally de-excited, any 
gas clump of density above 500 cm~'^ will induce coUi- 
sional de-excitation of the infrared cooling lines, boosting 
the [O III] A5007 line. Another problem is that in many 
cases photoionized gas at temperatures below ~ 5000 K is 
thermally unstable. Thus the temperatures computed in 
photoionization models do not reflect the whole range of 
temperatures that can be found in real objects. For these 
reasons the predicted intensities of the models with Z/Zq 
> 1 (corresponding to 0/H = 8.51 x 10"'' as stated above) 
are less accurate than those of lower metaUicity models. 
Another notorious problem of high metaUicity models is 
that the intensities of optical forbidden lines, especially 
from the low ionization species, are very dependent on 
the amount of metals trapped in dust grains (Henry 1993, 
Shields & Kennicutt 1995). Our reference grids of models 
assume no depletion, although we have constructed addi- 
tional models with only one tenth of the metals in the gas 
phase, see Sect. 4.3) 

4. Observations versus models 

We first analyze diagnostic diagrams involving the H/3 
equivalent width. The classical, pure emission line ratio 
diagnostic diagrams will be discussed subsequently. 

4-.1. Diagnostic diagrams with equivalent widths 

The observed samples of H ii galaxies, and giant H ii re- 
gions in general, are characterized by a very narrow range 
in ionization parameters. This was already noted before 
from a comparison of observed diagrams with model grids 
(McCall et al. 1985, Dopita & Evans 1986, Garcia- Vargas 
& Diaz 1994, Garcia- Vargas ct al. 1995, Stasinska & Lei- 
therer 1996, Bresohn et al. 1999, Dopita et al. 2000). We 
show in Fig. 12 the values of [O ni]/[0 ii] as a function of 
EW{}if3) for models with = 1, 10^, 10^ and 10^ M© 
and metaUicity 0.2 Zq (corresponding to a range of a fac- 
tor 10 in U between each curve) superimposed on the Ter- 
levich data. For EW{}i(3) > 100 A, the observed disper- 
sion in [O iii]/[0 ii] suggests a range in ionization param- 
eters of about a factor 10. Most of the data are bracketed 
by the series with M*= 10^ and 10^ Mq (corresponding 
roughly [/ ~ 1 - 2 x lO'^ and to [/ ~ 1 - 2 x lO'^, for the 



first 5 Myr). At lower equivalent widths, the observations 
fall outside the predicted values, apparently indicating a 
much higher U. As we will show below, this would be a 
premature conclusion. 

If H II galaxies are powered by coeval star clusters, one 
expects the [O iii] A5007/H/3 ratio to decrease as time goes 
by and EW(Hf3) decreases. Panels a in Figs. 6-11 show the 
predictions for [O ill] A5007/H/3 as a function ofEW{lip). 
All the sequences show a steep decline of [O iii] A5007/H/3 
at an EW (H/3) around 20 A. irrespective of metaUicity. A 
similar effect was predicted by the models of SL96. This 
is due to the average effective temperature of the exciting 
stars (as roughly expressed by Qhc"/Qh) dropping much 
more rapidly than the total number of ionizing photons, 
once the most massive stars have disappeared. The rise in 
the curves in Figs. 6-11 prior to the final steep decline for 
the Z = 2 Zq and Zq models is due to the appearance 
of Wolf-Rayet stars. The sequences of models with M^p = 
30 Mq show a strong accumulation of points at early ages 
since it takes stars with masses below 30 Mq more than 
4 Myr to evolve off the main-sequence. The drop in [O III] 
A5007/H/3 at EW{Rl3) smaller than 20 A is of course the 
same as for the M^p = 120 Mq sequences. The sample of 
H II galaxies in the SL96 study had EW{B.(3) too large to 
show this drop, as a consequence of the requirement on 
[O III] A4363, as stated in the introduction. Since no such 
restriction was applied here, the values of EW{lif3) go 
down to a few Angstroms. The data in each sample have a 
tendency to a decrease of [O ill] A5007/H/3 with decreasing 
EW{li0). This decrease, however, is rather mild and the 
scatter in [O ill] A5007/H/3 increases for lower EW(H(3). 
One might think that extending the star formation over a 
longer period would help to reconcile the predicted trends 
with the observations. However, Fig. 9 suggests that even 
considering a constant star formation during 10 Myr is 
not sufficient to produce the large [O iii]/H/? observed at 
EW{R(3) of 20 A or smaUer. 

A similar effect is seen in the He I A5876/H/3 versus 
EWiUP) diagrams (panels b). In contrast to [O iii]/H/3 
which, for a given stellar radiation field depends on U and 
Z, He I A5876/H/3 is independent of the ionization param- 
eter and depends only weakly on metaUicity (in the sense 
that at larger metallicities a lower electron temperature 
increases He i A5876/H/3). Instantaneous starburst mod- 
els predict a sharp drop in He i A5876/H/3 at EW{Ill3) 
around 20 A. This is not seen in any of the observed sam- 
ples. The observed He i A5876/H/3 ratio in the Terlevich 
sample remains remarkably constant around a value of 
0.1. The slight rise and the increasing dispersion at low 
EW{B.0) could be due to observational errors, He i A5876 
being a weak line. Unfortunately, this line is measured in 
only 98 out of 305 objects, and one might argue that the 
remaining objects have a much lower He i A5876/H/3 ratio. 
All the objects in the Izotov sample have the He i A5876 
intensity measured with very good accuracy, and all ex- 
cept one are of the order of 0.1. Although the statistics are 
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very poor, since only 10 objects in the Izotov sample have 
EW(RP) of the order of 20 A or smaller, there seems to be 
no drop of He i A5876/H/3 at this EWiUP). The Popescu 
& Hopp sample shows the same behavior as the Terlevich 
sample. 

The presence of an older underlying stellar component 
in H II galaxies, as supported by many studies based on 
infrared colors and stellar features, changes the interpre- 
tation of the diagrams. It would decrease the H/3 equiv- 
alent width of the young, ionizing population ( "dilution" 
effect), and also lead to an underestimate of the nebu- 
lar H/3 luminosity due to underlying stellar absorption. 
Qualitatively, these effects work in the direction required 
to reconciliate models and observations. Such an expla- 
nation is suggested by a recent study of Raimann et al. 
(2000a). These authors grouped 185 spectra from the Ter- 
levich et al. (1991) paper into 19 templates in order to 
increase the signal-to-noise and to allow a population syn- 
thesis analysis from the absorption features. They found 
that H II galaxies are age-composite stellar systems, and 
that they can contain a significant population of stars with 
ages larger than 10 Myr. In Fig. 13 we plotted the same 
quantities as in Figs. 1-3^ for the spectral groups defined 
by Raiman et al. (2000a) (excluding the 3 groups corre- 
sponding to Seyferts which are not discussed in the present 
paper). Open circles represent the averaged observed spec- 
tra. Each open circle is linked by a straight line with a 
filled circle representing the data after correction for an 
underlying stellar population (line and continuum) and 
for internal reddening. The same trends with EW{llf3) 
are seen for the filled circles as for the data shown in 
Figs. 1-3, but with larger slopes. The drop in [O iii]/H/3 
starts at much higher values of EW{}i(3). This is much 
more compatible with our models than the uncorrected 
data. The steep decline starts immediately, which seems 
to be in better agreement with our models with Mup= 30 
M0(Fig. 10) than with 120 Mq. However, the distribution 
of the open circles in the [O iii]/H/3 versus EW(RP) dia- 
gram of Fig. 13 is quite different from the distribution of 
the individual galaxies from the Terlevich sample (Fig. 1); 
therefore this averaged sample is unsuitable for describ- 
ing the general properties of H ii galaxies. The important 
point is that the work of Raimann et al. (2000a, b) demon- 
strates that signatures of older underlying stellar popula- 
tions are present in the spectra of H ii galaxies and can 
explain the observed trend in the uncorrected [O iii]/H/3 
versus EW{B.(3) diagrams of Figs. 1- 3. The multiwave- 
length analysis of 17 H li galaxies by Mas-Hesse & Kunth 
(1999) reached a similar conclusion. Using indicators of 
young stellar populations (the UV continuum, the stellar 
-I- interstellar W{Si w)/W{C iv) ratio, the strength of 
the Wolf-Rayet bump), these authors conclude that the 
majority of these galaxies have experienced recent nearly 
instantaneous star formation dominating the ultraviolet 



except He I A5876/H/3 which is not available 



light. Yet their synthetic spectra are often below the ob- 
served visible continuum, indicating the presence of an 
older stellar population. 

Other causes may contribute to the behavior of 
[O iii]/H/3 and He i A5876/H/3 as a function of EWiUP) 
as well. For example, H ii regions may become gradually 
density bounded or the covering factor may decrease with 
time. A lower covering factor decreases the H/3 luminosity 
and thus EW(Hf3). Allowing for leakage of ionizing pho- 
tons from the outskirts of the emitting region maintains 
[O III] /H/3 and He i A5876/H/3 at a high level. This is an at- 
tractive possibility, since as nebulae expand, they become 
less and less opaque to ionizing photons (for a homoge- 
neous gas sphere, the mass that can be ionized changes 
as Qn/n). Alternatively, it is possible that most of the 
giant H ii regions in the samples are density bounded or 
have a covering factor smaller than 1, regardless of their 
evolutionary state. Another possibility is that the spectro- 
scopic data encompass only part of the nebula (for sim- 
plicity, we will refer to these explanations as the "aperture 
effect" ) . Such "aperture effects" are likely to occur in our 
H II galaxy samples. A detailed study of the H ii galaxy 
I Zw 18 indicates that the brightest H ii region, which 
dominates the spectrum, is density bounded, at least in 
some directions (Stasihska & Schaerer 1999). This could 
be a general property since diffuse or filamentary ionized 
gas has been observed at large distances in a number of 
H II galaxies (e.g., Hoopes et al. 1996; Martin 1997; Hunter 
& Gallagher, 1997), implying that the data in this paper 
may be missing part of the distant low-ionization gas. 

High values of [O iii]/H/3 and He i A5876/H/3 at low 
EW{}if3) could also result if the extinction of the ionized 
gas is systematically higher than that of the stellar light, 
as found in several studies of H ii galaxies and starbursts 
(e.g. Fanelh et al. 1988, Calzetti 1997, Mas-Hesse & Kunth 
1999, Schaerer et al. 2000). 

To summarize, simple photoionization models of H ii 
galaxies as pure evolving starbursts cannot entirely ac- 
count for the observed diagrams involving EW(R(3) (pan- 
els a - d) of Fig 1-3. One or more of these effects are likely 
explanations: (i) an older stellar population; (ii) starlight 
systematically suffering a smaller extinction than ionized 
gas; or (iii) an aperture effect where the observed spectra 
correspond to density bounded nebulae or nebulae with a 
covering factor smaller than 1. These interpretations ex- 
plain in a natural way the absence of H ii galaxies with 
EW{ll/3) larger than 400 A in objective prism surveys, 
whereas starburst models for ionization bounded H li re- 
gions predict values of up to 700 A (see also discussion 
in SL96). This also implies that such surveys are biased 
against galaxies where the latest starbursting episode is 
older than about 5 Myr, as such objects would tend to 
have too weak H/3 and Ha emission lines to be detected. 

How can these explanations be reconciled with the in- 
crease of [O i]/H/3 as EW{lip) decreases (panels c in Figs. 
1-3)? This trend is also present in the corrected template 
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spectra from Raimann et al. (Fig. 13), but the number 
of points is small and confirmation on a larger sample is 
required. One could argue that the old underlying pop- 
ulation is stronger in higher metallicity objects, so that 
the increase in [O i]/H/3 with decreasing EW{HP) could 
be attributed to a metallicity effect. We can test this idea 
by plotting [O i]/([0 n] + [0 ill]) as a function of 0/H for 
the 60 objects with measured [O ill] A4363 of the Izotov 
sample. This is done in panel b of Fig. 14. Panel a of this 
figure shows the same objects in the [O i]/([0 11] + [O ill]) 
vs. EW(H/3) plane. Clearly, the metallicity effect is not 
a good explanation. Also, the behavior of [O l] cannot be 
explained by the "aperture effect" discussed above. Reduc- 
ing the covering factor while keeping the nebulae ioniza- 
tion bounded does not change the [O ii]/H/3 and [O i]/H/3 
ratios. Making the nebulae density bounded reduces the 
ratios. Note that the observed behavior of [O i]/H/3 with 
EW{B.p) argues in favor of EW(Rf3) remaining a statisti- 
cal age indicator (although the age cannot be obtained by 
a direct comparison with photoionization models for pure 
starbursts). As pointed out by by Stasihska & Schaerer 
(1999) in the case of I Zw 18, strong [O i] emission is 
easily produced by photoionization models in dense fila- 
ments which may contribute only marginally to the total 
emission in the other lines. Possibly, giant H 11 regions 
are globally density bounded but contain high density fil- 
aments. Such filaments could be produced by shock-wave 
compression or instabilities, produced by stellar winds and 
supernovae whose effects are expected to increase with 
time. 

4-2. On the strong line method to derive oxygen 
abundances 

An important consequence of these observational selec- 
tion effects in H 11 galaxy samples is that the age of the 
ionizing population has little influence on the observed 
([O iii]-f [O ii])/H/3 ratio (see panel e in Fig. 6-9). Com- 
bined with the narrow range of U in giant H 11 galax- 
ies (cf. panels a of Figs. 1-3 and Figs. 6 - 11), this 
explains a posteriori why strong line methods based on 
([O iii]4-[0 ii])/H/3 and methods utilizing the electron 
temperature lead to similar oxygen abundance estimates 
in these objects (Pilyugin 2000). A priori one would have 
expected large differences induced by the temporal varia- 
tion of the ionizing radiation field during the evolution of 
a star burst. 

Fig. 15 shows how the [O iii]/[0 11] vs 
([O iii]-t-[0 ii])/H/3 diagram (introduced by McCaugh 
1991, 1994 to derive 0/H and U at the same time) 
changes over a period of 4 Myr. Continuous lines join 
models with equal metallicities (the thickest lines cor- 
respond to the most metal rich cases), dotted lines join 
models with equal A/*. Panels a, b, c correspond to 
starburst ages of .01, 2.01 and 4.01 Myr, respectively. 
The objects from the Terlevich sample are overlaid in 



the three panels. We see that the theoretical diagrams 
at metallicities lower than half solar evolve little over 
this period, demonstrating the applicability of the strong 
line method at these metallicities. In practice, it is 
however preferable to calibrate the method not on models 
but empirically with a sample of high signal-to-noise 
observations of H 11 galaxies, as done by Pilyugin (2000). 

Interestingly, as already pointed out by Stasihska 
(1999), rapid evolution of the McGaugh (1991) diagram is, 
however, predicted at metallicities ^ solar, due to impor- 
tant changes in the hardness of the stellar radiation field at 
such metallicities, and due to the strong impact on optical 
line ratios. This behavior of the stellar radiation field at 
large metallicity is due to the fact that the overall main se- 
quence is shifted toward lower temperatures, and thus the 
temperature of the stars falls below Teg ^ 40 kK earlier, 
as the hardness measured by Qhb/Qh decreases rapidly. 
This explains the more rapid decrease of the hardness of 
the ionizing flux of stellar populations at high metallicity 
(cf. Schaerer & Vacca 1998, Leitherer et al. 1999)|]. There- 
fore at metallicities around solar or above (where electron 
temperature based methods cannot be used), the results 
from strong line methods to derive the oxygen abundance 
are likely to be rather inaccurate. 

We also emphasize that strong line methods are statis- 
tical, as compared to electron temperature methods. They 
should be used with caution in cases when a systematic 
variation in the starburst properties or the gas density 
distribution is expected. This applies to studies of envi- 
ronmental effects on metallicities, e.g., such as those of 
giant H 11 regions seen in tidal tails of galaxies (Due & 
Mirabel 1998) or H 11 galaxies in the core of galaxy clus- 
ters (Vilchez 1995) as compared to isolated H 11 galaxies 

4-3. Pure emission line ratio diagnostic diagrams 

Next we consider classical diagnostic diagrams which ex- 
clusively rely on line ratios. To a first approximation, such 
diagrams depend only on the population of the most mas- 
sive stars which produce the ionization. They have been 
extensively studied for giant H 11 regions in spiral galax- 
ies and, as mentioned in Section 2.4, it has been found 
that [O iii]/H/3 versus [O ii]/H/3, or [O iii]/H/3 versus 
[O 11] /[O III] define a very narrow sequence, called the 
H II region sequence. This sequence was analyzed with 
the help of photoionization models and has been inter- 
preted as a sequence in metallicity with effective temper- 
ature/ionization parameter varying in line with metallicity 
(McCall et al. 1985, Dopita & Evans 1986, Dopita et al. 
2000). 

^ Although in general the hardness decreases with age, the 
temporary increase in the high metallicity models at 4 Myr is 
due to the appearance of WR stars at ^ 3 Myr (cf. Fig. ^). 
The reality of this feature is, however, questionable (cf. Section 
3.1). 
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Since H II galaxy samples are believed to contain 
mostly low-metallicity objects, it is not surprising that 
the H II galaxy sequences shown in Figs. 1 - 3 a are less 
populated at the low [O ili]/H/3 end than the McCall et al. 
(1985) or van Zee et al. (1998) sequences, and that they 
extend to high [O ili]/H/3. 

Although some of the observed trends are met by the 
predictions, it is, nevertheless, apparent from Figs. 6 - 
11 that our models do not reproduce the sequence very 
well. Compared to other classical diagnostic diagrams, the 
[O lii]/H/3 vs [O ii]/H/3 diagram has the advantage of be- 
ing independent of the abundances ratios of 0/N or 0/S. 
Therefore we will focus our discussion on this diagram, but 
the conclusions are similar for the [O iii]/H/3 vs [S ii]/H/3 
diagram. Diagrams involving [N ll] are discussed in Sect. 

Excluding the 2 and Zq sequences (models repre- 
sented by squares and circles), which are not representa- 
tive of the bulk of the objects in our samples, the param- 
eter space occupied by the models in the diagram does 
not exactly match the observations. First, our model se- 
quences predict an almost vertical drop in the [O iii]/H/3 
vs. [O ii]/H/3 diagram, and the maximum [O ii]/H/3 is 
almost independent of the ionization parameter. This is 
due to oxygen mainly being in the form of O"*" when the 
most massive stars have disappeared. Since we have ar- 
gued above that our H ii galaxies sample should contain 
few, if any, galaxies with the most recent starburst older 
than 5 Myr, our observational diagrams should be com- 
pared to the model sequences until 5 Myr only (roughly 
corresponding to EW(HP) ^ 50 A in Figs. 6 - 11). 

Even with that restriction, our photoionization models 
are unable to reproduce [O ll]/H/3 ratios larger than 4 and 
the kink at [O iii]/H/3 > 4 and [O ii]/H/3 > 2. Moreover, 
the problem may become even worse if the model spectra 
used for the Wolf-Rayet phase were too hard (cf. Section 
0). 

We have computed supplementary photoionization 
models especially for this purpose (including depletion of 
metals on grains or heating by X-rays), but without suc- 
cess. The disagreement becomes even larger if we assume 
that H II galaxies are density bounded, as argued above. 
Whether additional heating sources (shocks, conduction, 
turbulence) may solve the problem is an open question. 

Dopita et al. (2000) compared an observed sample of 
giant regions in spiral galaxies with photoionization mod- 
els and found that their models could not reproduce the 
high [O i]/H/3 and [S ii]/H/3 ratios seen in some objects, 
invoking shocks as an explanation. Our models use a dif- 
ferent prescription for the density structure than those of 
Dopita et al. and differ also in other respects, such as nu- 
merical aspects and the N/0 vs. 0/H prescription. Still, 
we find that the problems exists not only for [O i]/H/3 and 
[S ii]/H/3 but also for [O ii]/H/3. Of course, the observed 
[O ii]/H/3 ratio is largely affected by the dereddening pro- 
cedure, while this is much less the case for the other two 



ratios. Since strong line diagnostics heavily rely on [O ii], 
and since this line is becoming increasingly important for 
the study of galaxies at high redshift, it would be useful to 
better constrain the problem by a detailed observational 
and theoretical study of those few H ii galaxies with a firm 
indication of shock excitation. 

McCall et al. (1985) argued that giant H ii regions 
must be mostly ionization bounded, otherwise one would 
not observe such a tight sequence in emission line ratio 
diagrams. This is at variance with the fact that we now 
know that a fair portion of Lyman continuum radiation 
is leaking out of H ii regions (Martin 1997; Stasihska & 
Schaerer 1999, Beckman et al. 2000). However, if the den- 
sity structure of giant H ii regions were driven by the me- 
chanical action of winds and supernovae explosions from 
their embedded stellar populations, one could understand 
why there is so little dispersion among giant H ii regions in 
such diagrams, even if they are partially density bounded, 
since the driving parameter would boil down to the stellar 
population itself. In this respect, hydrodynamic modelling 
of the ionization structure of giant H II regions (the follow- 
up of the single star H ii region modelling performed by 
Rodriguez-Caspar & Tenorio-Tagle 1998) would be ex- 
tremely important. 

4.4. The N/0 ratio 

Among others, Kunth & Sargent (1986), Pagel et al. (1986, 
1992), Olofsson (1995b) suggested that gas might be chem- 
ically enriched in, e.g., helium and nitrogen from the winds 
of Wolf-Rayet stars during the evolution of giant H 11 re- 
gions. Traces of such a local pollution have been found 
in the irregular galaxy NCC 5253 from imaging spec- 
troscopy (Walsh & Roy 1987, 1989, Kobulnicky et al. 
1997), but other attempts to detect local N/0 enhance- 
ment in Wolf-Rayet galaxies have failed (Kobulnicky & 
Skillman 1996, 1998, Kobulnicky 1999). Moreover, Izotov 
et al. (1997) and Kobulnicky & Skillman (1996) analyzed 
observational data on H 11 galaxies in the literature and 
found that galaxies with strong Wolf-Rayet features in 
their integrated spectra exhibited the same N/0 ratios as 
the remaining galaxies at identical 0/H ratios, contrary 
to the claim of Pagel et al. (1986, 1992). 

In the light of these previous studies, the [N 11] /[O 11] 
vs EW{llf3) diagram is revealing. There is a definite trend 
of [N ii]/[0 11] strongly increasing as EW{H0) decreases 
in all our H 11 galaxy samples. If we take [N 11] /[O 11] as 
an indicator of N/0, it is tempting to attribute this trend 
to an increase in N/0 as the ionizing starburst gets older. 
However, several biases must be examined first. Could the 
observed relation be the result of a selection effect? While 
the absence of low [N 11] /[O 11] ratios at small equiva- 
lent widths could be attributed to selection effects against 
the weakest lines, at large H/3 equivalent widths, values of 
[N 11] /[O 11] larger than ~ 0.2 - 2 should be observed if such 
objects exist. Moreover, there is no selection effect in the 
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Izotov sample since all the objects appear in Fig. 2f, and 
the relation is seen there as well. The trend is less distinct 
than in the other samples, probably because the sample is 
weighted towards the most metal poor galaxies. Finally, 
as mentioned in Sect. 2.4, if a complete subsample is ex- 
tracted from the Terlevich sample, the observed relation 
between [N ii]/[0 ii] and EW{}ip) remains. Therefore, the 
observed relation cannot be attributed to a simple selec- 
tion effect. 

Turning to an interpretation, it is important bearing in 
mind that [N ii]/[0 ii] does not only depend on N/0 but 
also on the electron temperature. It is larger at lower elec- 
tron temperatures and therefore at higher O /H for a given 
N/0. Second, at metallicities above solar, the [N ii] and 
[O ii] lines are produced by recombination rather than by 
collisional excitation (if the N++/N+ and 0"'""'"/0"'" ratios 
in the nebula are not close to zero). Our models (panels 
f in Figs. 6 - 11) do account for that. We have already 
argued that models with ages ^ 5 Myr may not be rel- 
evant to our samples of H ii galaxies, but that EW{B.p) 
is an indicator of the age of the ionizing starburst (in the 
sense that smaller EW(H(3) correspond to larger ages). 
We have also argued that one should perhaps consider 
density bounded rather than ionization bounded models. 
However, the [N ii]/[0 ii] ratio would remain imchanged in 
a density bounded model. The large values of [N ii]/[0 ii] 
(^ 0.3) seen in some of the H ii regions of our samples 
therefore indicate a N/0 ratio larger than the largest pre- 
dicted by our models (0.1 at Z = 2 Z©). 

If the metallicities (i.e., O/H ratios) in our samples 
have identical distributions in each age bin. the observed 
trend could be explained if the two following conditions 
hold at the same time: (i) the relation between N/0 
and O/H is steeper than the relation based on McGaugh 
(1991) adopted in our models and (ii) the most metal rich 
objects are those whose EW{li(3) values are the most af- 
fected by the contribution of an old stellar population. 
The first hypothesis might indeed be true. There are in- 
direct indications, from studies of giant H II regions in 
spiral galaxies, that at high metallicities N/0 increases 
more rapidly than N/0 a (0/H)°-^ adopted in our mod- 
els (cf. van Zee et al. 1998, Henry et al. 2000). The second 
hypothesis appears quite reasonable. Earlier generations 
of stars are necessary to produce the bulk of presently 
observed O/H in more metal rich H ii galaxies, and these 
earlier generations still contribute to the continuum. Qual- 
itatively, the results of Raimann et al. (2000b) support 
this picture. We note however that the anticorrelation be- 
tween [N ii]/[0 ii] and EW(Hl3) subsists in the Raimann 
et al. templates even after correction for the old stellar 
population (Fig. 13f). If this effect is real, it calls for an 
additional explanation. 

N/0 might be increasing with time. Such an ex- 
planation would be more compatible with the observed 
[N ii]/[0 ii] versus ([O ii] + [0 iii])/H/3 diagram in Figs. 1 - 
3 i, in which the scatter in ([O ii]-|-[0 iii])/H/3 at a given 



[N il]/[0 ii] is extremely small. Alodcls with ages smaller 
than 5 Myr (panels i in Figs 6-11) show about the same 
scatter, but in order to compare with observations a con- 
volution of the theoretical scatter with observational er- 
rors, and a correction for reddening and underlying H/3 
absorption is first required. With N/0 increasing with 
time, the theoretical scatter at a given [N ii]/[0 ll] be- 
comes smaller. The hypothesis of N/0 increasing with 
time can in principle be checked directly by plotting the 
N/0 determined by electron temperature based methods 
as a function of EW{}i(3). We performed this test. No 
obvious trend is seen, but none of the objects with large 
observed [N ii]/[0 ii] in our samples has [O III] A4363 
measured, so this test is not conclusive. 

Deep spectroscopy for a detailed analysis of the stellar 
populations in the objects showing high [N ii]/[0 ii] ratios 
should be obtained. High resolution emission line imaging 
and tailored photoionization modeling should also be un- 
dertaken, in order to put the strongest possible constraints 
on O/H and N/0 in these objects. 

5. Conclusions 

We have reanalyzed the emission line properties of three 

large samples of H II galaxies taken from the literature 
(Terlevich et al. 1991; Izotov and collaborators 1994-2000; 
Popescu & Hopp 2000) with the aim of studying the tem- 
poral evolution of these objects during the lifetime of the 
ionizing stars (^ 10 Myr). We have constructed a series of 
diagrams using observed emission line ratios and equiva- 
lent widths, and found significant trends. 

We interpreted these trends with photoionization mod- 
els for integrated stellar populations. This study extends 
the work of Stasinska & Lcitherer (1996), by including ob- 
jects with no direct metallicity determination and using 
updated evolutionary tracks and atmospheric models. 

The changes brought about by the inclusion of the 
latest non-rotating Geneva stellar evolution models and 
the CoStar non-LTE line blanketed atmosphere models 
including stellar winds for O stars have been summarized 
in Sect. 3.1. (cf. Fig. 5). These models have been suc- 
cessfully compared to observations in a large variety of 
studies related to massive stars (sec e.g. Maedcr & Meynet 
1994, Stasinska & Schaerer 1997, Oey et al. 2000, Schaerer 
2000). Despite this, the adopted models have also some 
known weaknesses, e.g. the neglect of stellar rotation lead- 
ing to modifications of the stellar tracks (cf. Maeder & 
Meynet 2000), a possible overestimate of the hardness of 
the ionizing flux above ^ 40 cV (Oey et al. 2000), or the 
lack of line blanketing in WR atmospheres (see Sect. 3.1). 
Essentially, the main dependence of the photoionization 
models on the stellar input physics enters through 1) the 
evolutionary timescales and the temperatures of main se- 
quence O stars, and 2) the hardness of the He°/H ioniz- 
ing spectrum. It is important to note that, while detailed 
outputs from photoionization models are sensitive to the 
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exact stellar input physics, our main conclusions summa- 
rized below are quite robust. 

The underlying principle of the comparison between 
models and observations is that the mctallicity distribu- 
tion is age-independent for sufficiently large and homoge- 
neous samples, despite the fact that the metallicities of 
individual objects may be a priori unknown. 

We have shown that H il galaxies selected by objective- 
prism surveys do not correspond to the simplistic view of 
instantaneous starbursts surrounded by constant density, 
ionization bounded H il regions. The observed relations 
between emission line ratios and H/3 equivalent width can 
be understood if most H ii galaxies contain older stellar 
populations contributing, sometimes rather significantly, 
to the observed optical continuum. This is in line with 
the recent findings of Raimann et al. (2000ab) who, from 
a spectroscopic analysis on the stellar features in very 
high signal-to- noise templates of H ii galaxies, showed that 
populations as old as 100 Myr and up to a few Gyr are 
detectable in the spectra and significantly afi^ect the ob- 
served H/? equivalent widths. This is also consistent with 
the model of a long lasting low level of star formation in 
I Zw 18, as suggested by Legrand (1999). Two additional 
mechanisms may play a role. First, stars seem to suffer 
a smaller dust obscuration in the visible light than the 
ionized gas. Second, some Lyman contimmm radiation is 
probably leaking out of most of the nebulae encompassed 
by the observing slits. 

The combination of these effects reduces EW{ll/3) 
with respect to the value predicted for an instantaneous 
starburst surrounded by an ionization bounded nebula. 
Therefore objective-prism selected samples of H ii galaxies 
are unlikely to contain significant numbers of starbursts 
older than about 5 Myr. Older stages in the evolution 
of starbursts must be selected from photometric surveys 
based on broad- or narrow-band colors. 

An interesting consequence of this selection effect is 
that the strong line methods for deriving oxygen abun- 
dances work rather well in metal poor H II galaxies be- 
cause there is no large mean effective temperature spread. 
Under these conditions, a strong lino method which is caH- 
brated on objects with electron temperature based oxygen 
abundances determinations like that proposed by Pilyugin 
(2000), is the most valid approach to derive oxygen abun- 
dances in low signal-to-noise spectra of H ii galaxies. One 
must, however, keep in mind the statistical nature of this 
method and the fact that the unknown stellar absorption 
in H/3 provides a further source of uncertainty. 

We have shown that [O i]/H/3 uniformly increases with 
decreasing EW{liP). This behavior suggests that this line 
ratio results from dynamical effects that shape the nebula 
and whose importance increases with time. These dynam- 
ical effects could also be responsible for the small range in 
ionization parameters that account for the observed emis- 
sion line trends. 



The classical emission line ratio diagnostic diagrams 
such as [O III] /H/3 vs [O ii]/H/3 imply a sequence in oxy- 
gen abundance and ionization parameter. This suggestion 
was made earlier on the basis of single-star photoionization 
models, but a physical interpretation remains elusive. Ex- 
tra heating sources in addition to the Lyman continuum 
radiation from massive stars seem necessary in order to 
explain the largest observed [O ii]/H/3 ratios. 

The [N ii]/[0 ii] ratio is shown to increase as EW(Rl3) 
decreases. The observed trend is even stronger when con- 
sidering the sample of high signal- to- noise H ii galaxy tem- 
plates of Raimann et al. (2000ab) after correction for the 
old stellar population. We conclude that cither the relation 
between N/0 and 0/H is steeper than that adopted in our 
models (N/0 oc 0/H°-^) and the underlying stellar popu- 
lation is stronger at higher metallicities, or the N/0 ratio 
increases with time on a time scale of ~ 5 Myr. This last 
option would support a scenario of self-pollution of giant 
H II regions by nitrogen produced in situ. High signal-to- 
noise observations to directly uncover the old stellar pop- 
ulations and the H/3 contamination in individual objects, 
especially those with the highest [N ll]/[0 ll] ratios, are 
required to settle this important issue. 

The possible dependence of the above findings on the 
various main unknowns (underlying populations, escape 
of Lyman continuum photons, non stellar heating sources, 
etc.) affecting our analysis have been discussed in the re- 
spective sections. Luckily, not all the problems raised in 
this study are interconnected, and they may be attacked 
from various angles. E.g. our claim of underlying popula- 
tions being responsible for an effective "age bias" against 
burst events with ages ^ 5 Myr and the puzzling behav- 
ior of [N ii]/[0 ii] can be tested by additional and more 
detailed stellar population studies of H ii galaxies along 
the lines of Raimann et al. (2000ab), using high signal- 
to-noise spectroscopy. A purely "stellar" solution seems 
now clearly excluded for the problem of [O iii]/Hi'3 vs 
[O ii]/H/3 as well as [S ii]/H/3 (Sect. 4.3 and Stasinska 
& Schaerer 1999). Hydrodynamical models of the interac- 
tion of the interstellar medium with clusters of hot stars 
(in the vein of Canto et al. 2000 or Franco et al. 2000) 
combined with photoionization calculations, should allow 
one to propose quantitative solutions. Such an approach 
may at the same time explain the behavior of [O i]/H/3 
as a function of EW{H.(3) and provide some quantitative 
estimate of the leakage of ionizing photons from giant H ii 
regions. Multi-wavelength high spectral resolution obser- 
vations should provide important constraints to the mod- 
els in this respect. The dust obscuration issue is perhaps 
trickier, since dust effects are extremely dependent on ge- 
ometry. Here, systematic high resolution imaging of giant 
H II regions, such as performed by Calzetti et al. (1997) or 
Johnson et al. (2000), should allow one to clarify the dust 
location and its effect on integrated H ii galaxies spec- 
tra. The insight gained from this study and the proposed 
directions for further theoretical and observational inves- 
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tigations will considerably increase our understanding of 
the physics of H ii regions and the complex interactions 
between its stellar and interstellar components. 
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Fig. 4. Emission line ratio diagrams for the sample of giant H II regions in spiral galaxies from Mc Call et al. (1985). 




Fig. 5. Comparison of the predicted hardness of the ionizing spectra from various evolutionary synthesis models for 
instantaneous bursts, Salpeter IMF, and different metallicities. Illustrated are the ratios of the He" and He+ to H 
ionizing photon fluxes (Qhc/Qh and Qne+Z^H respectively), a) Schaerer & Vacca (1998, hereafter SV98) models 
(solid, long-dashed) versus Leitherer & Heckman (1995) models (dashed, dotted) at solar metallicity. b) SV98 models 
(solid, long-dashed) at 1/5 solar, Leitherer & Heckman (1995) models (dashed, dotted) at 1/4 solar, c) SV98 models 
(solid, long-dashed) versus Starburst99 models (dashed, dotted) at 1/5 solar, d) Same as c for 1/20 solar metallicity. 
Discussion in text. 
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Fig. 9. Same as Fig. 7, Afup = 120 Mq, M*= 10^ Mq, but for an extended burst (duration 10 Myr). Same symbols 
as in Fig. 7 but with time steps of 2 Myr. 
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Fig. 11. Sequences of models for instantaneous bursts (Afup = 120 M©, M+= 10^ M0). hollow sphere geometry, n 
= 200 cm~^, various metallicities. These models have the same average ionization parameter as the models of Fig. 7. 
Same symbols as in Fig. 5. 
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Fig. 13. Groups of H ii galaxies as defined by Raimann ct al. (2000a) from the Terlcvich ct al. (1991) catalogue. Open 
circles: the averaged spectra as used by Raiman et al. (2000a); filled circles: the same after correction for underlying 
stellar population and internal reddening. 
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Fig. 15. Evolution of the McGaugh (1991) diagram over a period of 4 Myr for our instantaneous burst models (M^p = 
120 Mq), full sphere geometry, n = 10 cm~^. Models with equal metallicities are linked by full lines (whose thickness 

is proportional to the metallicity). Models with equal are linked with dotted lines. Panels a, b, c correspond to 
starburst ages of .01, 2.01 and 4.01 Myr, respectively Superimposed are the objects from the Terlevich sample 



